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Abstract 

We describe the spatial structure of particles in the (one dimensional) two-species annihilation 
reaction A + B — » 0, where both species have a uniform drift in the same direction and like species 
have a hard core exclusion. For the case of equal initial concentration, at long times, there are 
three relevant length scales: the typical distance between similar (neighboring) particles, the 
typical distance between dissimilar (neighboring) particles, and the typical size of a cluster of one 
type of particles. These length scales are found to be generically different than that found for 
particles without a drift. 

PACS numbers: 5.40.+j, 82.20.Mj, 2.50.-r 

1 Introduction 

In the irreversible two species annihilation reaction A + B — ► 0, local fluctuations in the particle 
density are sufficiently important that (at least in low dimension) the mean field rate equations do 
not accurately predict the decay of the density when the initial concentrations of the two species are 
equal [1, 2, 3, 4, 5, 6, 7, 8] 

There is a standard picture for what does happen [1]: in a region of size L it typically takes a 
time on the order of Z 2 for all of the different particles to react, since they may need to diffuse across 
the entire region in order to annihilate. One species in this region will be completely eliminated and 
excess particles of the other species will remain; their number is the initial excess of either type A 
or type B particles, and its expected value is proportional to the square root of the volume. At time 
t the concentration c(t) thus behaves as 

c(t) ~ (V) 1/2 jh d ~ (V/ 2 ) 1/2 jt d i 2 = r d i\ (i) 

for dimension d < 4, i.e. t -1 / 4 in dimension one, which has been verified rigorously [4, 7]. 

The above result depends crucially on the fact that the particles undergo ordinary diffusion. 
Yet this might not be correct if there is a drift field, even if both species drift in the same direction. 
Clearly one expects to see different behavior if the two species drift in opposite directions: in addition 



to the diffusive length scale L± ~ i 1 / 2 there will be a drift length scale L| ~ ut; the simple heuristic 
reasoning that lead to eq. (1) in the drift-free case then yields [9, 10] 

c(t) ~ (^( rf+1 )/ 2 ) V2 / (yt {d+1 ^ 2 ) = v l/2 < -(d+D/4 (2) 

for dimension <i < 3, or (v/t) 1 ^ in one dimension. 

In a previous work [8], we claimed the more suprising result that even if the two species have 
a drift in the same direction, even if their relative velocity is zero, the result may still be different 
from the expected t~ d ' 4 . This is because after subtracting off the average motion, there may be 
a super- diffusive length scale, of order t 2 ' 3 in one dimension [11, 12]. Thus one might expect the 
concentration to go as 

C(t) -I 1 ' 2 /!,- [ t ^] V2 /i 2 / 3 = i-1/3 . ( 3) 

This occurs when the long time behavior of the system is governed not by a linear diffusion 
equation, but by a nonlinear equation, the noisy Burgers equation [13], 

— + p— = v— + — (4) 

dt dx dx" 1 dx ' 

where p is a rescaled density, v is the viscosity (representing the lattice spacing) and £ is a random 
noise term, e.g. white noise where the covariance is 

{£(x, t)$(x', t')) = 2v6{x - x')6{t - t') . (5) 

As an example of a system that falls into this universality class we choose the Asymmetric Simple 
Exclusion Process (ASEP). We thus consider 2 species of particles that individually undergo ASEP 
dynamics but annihilate when jumping on top of each other. 

Both the drift and the exclusion rule are necessary for Burgers equation (4) to be relevant; a drift 
without the exclusion rule has no effect on the asymptotic behavior of the annihilation dynamics [14]. 
To keep matters simple we consider infinite drift — particles move either to the right or not at all. We 
use periodic boundary conditions on the one-dimensional lattice. A short description of the dynamics 
follows: at each time step, a site is picked at random. If that site is occupied, its particle attempts 
a jump to its nearest neighbor on the right. If the target site is empty, the jump succeeds. If the 
target site is occupied by a particle of the same species, the jump fails. If the target site is occupied 
by a particle of the opposite species, both particles annihilate. 

The asymptotic behavior of the density in this system was studied in [8]. Here we provide further 
evidence that this system is in a different universality class than the system without a drift by 
examining the spatial organization of the reactants in this system. This type of analysis was done 
by Leyvraz and Redner [5, 6] for the drift-free model; our results are unambiguously distinct from 
theirs. 
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Figure 1: Log-log plot of the average same-species interparticle spacing, (lAA)t- The dotted line has 
a slope of 1/3. 



2 Distribution of interparticle distances 

As in [5, 6], our most important results deal with three distributions: 1) Pt{lAA = x ), the probability 
that nearest neighbor particles of the same type are at distance x at time t\ 2) Pt(lAB = x ), the 
probability that nearest neighbor particles of different type are at distance x; 3) Pt(L = x), the 
probability that a domain of one type of particle has length x. The length of a cluster L is defined as 
the distance between the first particles of opposite species on either side of the cluster. We performed 
numerical simulations to measure these quantities; our results are based on averaging 9 independent 
runs of a system of size 4,000,000; times ranged up to 512,000 time steps. The initial density of each 
species was 0.45. (In [8] this density seemed to give the fastest approach to the asymptotic regime.) 

We first consider Pt(lAA = x )- F° r & U times t, this appears to be a monotonically decreasing 
distribution, almost exponential in shape. The average same-species interparticle spacing increases 
with time, as is shown in figure 1; the average spacing increases as t 1 ' 3 . Recall that in the symmetric 
model the spacing grew as i 1 ' 4 . Here the interparticle spacing grows more rapidly; consistent with 
the more rapid decrease in density (i -1 ' 3 ) observed in [8]. It is interesting that the density took 
some time to reach the asymptotic regime, while the spacing seems to obey scaling even for very 
short times. 

An example of the distribution itself is plotted in figure 2; the time here is 8000 but the graphs 
for other values of t are similar. It appears to be the superposition of two exponential distributions — 
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Figure 2: Semi-log plot of the (unnormalized) probability distribution for the same-species 
interparticle spacing at time 8000, Pti^AA = x )- The data is smoothed by averaging the distribution 
over 5 consecutive points. The slope of the steeper straight line is exactly double that of the other. 



one with twice the decay rate of the other. The explanation for this is quite simple: when clusters 
coalesce, the new cluster that is formed is the sum of the lengths of the constituent clusters. Thus, 
neglecting correlations between neighboring clusters, the distribution is (approximately) proportional 
to 

exp(— x/{lAA)t) + reaction rate x exp(— x/{2lAA)t) , (6) 

since the clusters which come together can be of arbitrary size. This same type of distribution is 
seen in the symmetric case, except of course for the different behavior of the average spacing. 

The distribution for the different-species interparticle spacing, or equivalently, the distance 
between domains, has a different shape entirely Since A and B type particles annihilate if they come 
too close, there is an effective repulsion between particles of different species. Thus Pt(lAB = x) is 
zero at x = 0, increases to a peak near {lAA)t (apparently linearly), and then decays (apparently 
exponentially); however, there does not seem to be a nice fit to the simple form x exp(— x/Xab) as 
is seen in the symmetric case [6]. 

Because of the effective repulsion, one expects that 

{lAA)t < (lAB)t ■ (7) 
In fact we find that not only is the distance between unlike particles greater than that between like 
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Figure 3: Log-log plot of the average different-species interparticle spacing, (lAB)t- The dotted line 
has a slope of 3/8. 

particles, but the growth rate is greater as well. Surprisingly, it appears to be the same growth rate 
as in the drift-free case: namely (lAtl)t ~ t 3 ' 8 , as is seen in figure 3. 

While 3/8, the opposite species growth rate, is only slightly larger than 1/3, the same species 
rate, this difference in length scales does mean that our clusters are in fact really clusters: inter- 
cluster distances are much larger than intra-cluster distances (at least for t large); this validates the 
assumptions of [8] which are based on treating clusters as isolated from each other. The fact that the 
ratio of these lengths scales only as i 1 ' 24 is perhaps responsible for the long times needed to observe 
asymptotic behavior for the density. 

3 Particle distribution within a domain. 

The third relevant length scale is the length of the clusters themselves. Since small clusters are quickly 
annihilated, and large clusters take a long time to form, we expect a single-peaked distribution 
for Pt(L = x) that is zero at x = 0. In fact qualitatively the distribution looks very similar to 

The average value for L as a function of time is plotted in figure 4. As in the drift-free case, the 
slope of this line can be determined from the other length scales in the problem [6]. Basically the 
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Figure 4: Log-log plot of the average domain size, (L)t- The dotted line has a slope of 7/12. 



concentration evolves via 

dp 



(8) 



&t (L)At ' 

where At is the time it takes for a nearest-neighbor pair of opposite species to annihilate, and (L) is 
the typical cluster length. 

If both members of the nearest-neighbor pair is assumed to diffuse, 

^diffuse ~ {Iab) 2 ■ (9) 

Note that this assumption of diffusive behavior is suspect: if we have an A-B pair, the A particle 
will move to the right with average velocity 1 and variance oc t. The B particle will do the same 
only so long as it does not hit another B particle. If collisions of the B particle are relevant, and 
if the density immediately to the right of the B particle is pboundary^ the B particle will move with 
average velocity 1 — ^boundary an d variance oc t^' 3 [15]. Thus the A particle slowly catches up to the 
B particle in a time 

Adrift ~ (U.b)/[1 - (1 - ^boundary)] = Qab) /^boundary , (10) 

provided that this is shorter than the diffusive time interval. 

Plugging in our results for the behavior of (Iab), and assuming that the boundary density is the 
same as the bulk, we have 

A^diffuse ~ Adrift ~ t 17 ^ 4 , (11) 
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Figure 5: Scaled domain profile at various times. L is the length of the cluster and x is measured 
from the center of each cluster. 



and thus, since p ~ t 4//3 , 

(L) t ~ t 7 / 12 or i 5 / 8 (12) 

depending upon whether drift or diffusion is relevant. Which is correct? The data in figure 4 seems 
consistent with the diffusive picture, i 7//12 , which can only occur if ^boundary *C Pbulk- 

How does ^boundary behave? Is ^boundary *C Pbulk, s0 that the observations and calculations are 
consistent? One way to approach this is to examine the density profile within a domain: in figure 5 
we plot the probability that a particle in a cluster of length I is a distance x/L from the center 
of the cluster. Unlike the case with no drift [6], this probability does not seem to have a scaling 
form. At best one could say that near the center of the distribution one is approaching a scaling 
form. In particular the behavior near the boundaries seems to be getting sharper and sharper. This 
suggests that perhaps x/L is not the correct scaling variable for cluster boundaries, i.e. the boundary 
behavior might scale differently from the bulk. 

A first guess might be that at the boundary between clusters the A-A distance might in fact scale 
similarly to the A-B distance. To check this we re-plot the data in figure 5 with a new horizontal 
axis of a:/ (Li -5 / 24 ), with x measured from the left edge of the cluster; this scaling is chosen because 
we want to rescale relative to 

(l^^tV^t^t-^-iLht- 5 ^; (13) 
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Figure 6: Scaled domain profile at various times. The position x is measured from the left edge of 
each cluster. 



when we do this in figure 6 we see that we obtain excellent scaling behavior over a very large range 
of times. This implies that t 3 / 8 is the correct length scale on the left-hand edge of the cluster, so 
that ^boundary ~ t~ 3 / s , which when inserted into eq. (10) gives 

^diffuse ~ At dri f t ~ t 17 / 24 , (14) 

and thus the result that 

{L) t ~t 7 ' 12 (15) 

is consistent with both derivations. 

If one instead looks at the profile near the right-hand edge of the cluster, neither (L)t nor {lAB)t 
appears to be the appropriate length scale. Of course because of the asymmetry of the system, one 
should not necessarily expect that the two boundaries should behave similarly. Nevertheless, one 
might still expect there to be scaling; perhaps with a different relevant length: is there any type of 
scaling that is correct here? 

The answer appears to be yes: although the scaling is certainly not as good as on the left-hand 
boundary, rescaling by Lt~ 12 seems to be appropriate), as seen in figure 7, where we plot the profile 
against Lt~' 12 . (A naiive error analysis would place the rescaling length between Lt~ A and Lt~' 14 .) 
While not entirely convincing, particularly given that we have no theoretical explanation for this 
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Figure 7: Scaled domain profile at various times. The position x is measured from the right edge of 
each cluster. 
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result, it is not entirely surprising that there is some means of superposing the data from different 
times. So apparently there is a fourth relevant distance in this mo del, t 7 l l2 t- A2 = t A6 . 

4 Conclusion 

In the ordinary, drift-free model for annihilating particles, there are three distinct (related) length 
scales. When a drift is added, we find four length scales. 

The first and most important length scale is (lAA)t ~ t 1 ^ . Since the clusters turn out to be large 
compared to the cluster separation, it is the nearest-neighbor same species particle distance that 
determines the density: p ~ 1/(1 AA)t ~ i -1 ^ 3 - 

For the inter-cluster distance we find that (lAB)t ~ t 3 / 8 . This is the only length scale unchanged 
by the presence of the drift field. 

The cluster size scales as (L)t ~ t 7 ' 12 . If this seems unlikely, recall that this is growth resulting 
from a combination of the expansion of individual clusters and the merging of pairs of clusters, so 
there is no reason for it to be particularly slow. 

Lastly, we have the length scale t 46 seen in the right-hand-side of each cluster. Because of 
symmetry, this length scale does not appear at all in the drift free case; its existence, as well as the 
different values for the exponents of the other length scales, indicates very clearly that the addition 
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of a drift changes the universality class for the problem of annihilating particles. 
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